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Abstract 

The transmission of fermions of mass m and energy E through an electrostatic potential barrier 
of rectangular shape (i.e. supporting an infinite electric field), of height U > E + mc 2 - due to 
the many-body nature of the Dirac equation evidentiated by the Klein paradox - has been widely 
studied. We exploit here the analytical solution, given by Sauter for the linearly rising potential 
step, to show that the tunnelling rate through a more realistic trapezoidal barrier is exponentially 
depressed, as soon as the length of the regions supporting a finite electric field exceeds the Compton 
wavelenght of the particle - the latter circumstance being hardly escapable in most realistic cases. 
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Figure 1: Klein step 
I. INTRODUCTION AND MAIN RESULT 

We will consider a one- dimensional flux of free monoenergetic electrons of energy E and 
momentum 

p = y/E 2 — m 2 (1) 

(throughout the paper natural units h—c—1 will be used) hitting a repulsive rectangular 
potential step of height U greater than their kinetic energy E — m (Fig. [T|) - to which they 
are minimally coupled. In this case total reflection is unavoidable: in the region x > there 
is indeed just one bounded solution of the Schrodinger equation and the flux is therefore 
zero. As a consequence, in the region x < the reflected flux must equal the incident one. 

The use of the Schrodinger equation is legitimate as long as both U, E — m <C m and, in 
principle, the result must not necessarily hold for higher values of U. Indeed for 

U — m > E > m (Klein zone) (2) 

(the only zone we will be interested in throught the paper) a relativistic equation is more 
suited for the description of the situation and, for both the Klein-Gordon and the Dirac 
equation (to which our discussion will be limited), it happens that in the region x > the 
plane-wave (free propagation) solutions with opposite momenta ±g, 

q = y/(U - E) 2 -m 2 (3) 

are two. Therefore the possibility of a non trivial transmitted flux of the same order of mag- 
nitude as the incident one is re-opened. The first to point out such a seemingly paradoxical 
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Figure 2: Sauter step 

result was Klein [l] in the case of the then recently proposed Dirac equation. The result 
has been commented upon and used by several authors over the years: see e.g. ref. Q for a 
historical review and some references. For us it is important to mention that Klein's result 
was questioned by Sauter [2I who, following a suggestion by Bohr, showed - at least for the 
Dirac equation - that, when the sharp edge of the step is substituted by a more realistic one 
of width t (FIG. E]), the transmission coefficient turns out to be 

^Sauter step ^ ^—-Krrfil/U 

The asymptotic form exhibited in (Tj0) holds when both the particles are fast E — m ~ m 
and the slope of the step, i.e. the electric field, satisfies 

e£ = U/£<^ 2m/A C ompton = 2m 2 . (5) 

A discussion, in the above framework, of whether the Klein paradox is a 'real' one and - at 
least for the Dirac equation - Sauter's is the right way out, cannot help spelling out what 
particles do the asymptotic states describe. Indeed, while in the region x < the dispersion 



relation E = + a/p 2 + m 2 is unambiguous, in the region x > £, owing to ([2]) and ([3]), either 

E = U - vV + m 2 (6) 



or 



- E = -U + vV + m 2 . (7) 

Namely either particles have negative kinetic energies, or antiparticles (propagating back- 
ward in time) have to come into play: the first choice being untenable, the many-body 
nature of both relativistic equations can no longer be ignored. 
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Figure 3: Klein barrier 

Postponing the discussion of this point until Section II, we prefer instead to examine the 
Klein paradox in the context of a different 'Gedanken' experiment where it is possible, for 
a while, to "sweep the dust under the rug": can an appreciable fraction of monoenergetic 
electrons pass through a potential barrier of height U > 2m ? 

The situation is summarized in FIG.s[3]andHl The former is the well studied rectangular 
barrier of width 2L whose transmission coefficient is known both in the Klein-Gordon and 
Dirac case: 

T-irect bar /q\ 

± ~ l + f ± (E,U) sin 2 (2gL) {) 



(U 2 (U/2-E) 2 _ + = RQ 



f±(E,U) 



p2 g2 



(9) 

m 2 U 2 



p2 g2 



D 



p, q being given by ([TJ), (J3J) respectively. This is the schematization of 'infinite electric field': 
the width of the edges of the barrier is neglected with an ensuing electric field much higher 
then its relevant scale 2m 2 /e given in (jSJ). This limitation is instead removed in the case of 
the trapezium shaped potential of Fig. HI the extension of Sauter's cure to the barrier. It is 
our choice to preserve space inversion as a symmetry (it will be evident that our conclusion 
does not critically depend on this assumption) and we will nickname this potential as Sauter 
barrier. 

The advantage of barriers with respect to steps is that in both the regions where {7 = 0, 
one may choose asymptotic states describing particles and ignore 'what is going on under 
the barrier' (we mean — L < x < L), i.e. whether either particles or antiparticles and/or 
couples are freely propagating. 
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Figure 4: Sauter barrier 

This aspect has been examined in Q] for the step, in [5| for the barrier, by adopting a 
space-time description of the scattering process instead of the stationary state picture we 
are sticking to here. 

n 

Also the authors of [5] claim a 'barrier paradox'. How does the paradox show up in 
the case of the rectangular barrier? Firstly it will be noted in ([8]) that for the values of 
E such that 2qL = nir, n integer, there is total transmission (the so-called transmission 
resonances), both in the Klein-Gordon and the Dirac case, no matter how high the values 
of U and L may be. This is in blatant contrast with the Schrodinger case mentioned in the 
beginning. A second crucial difference is in the role played by the barrier width parameter 
L. While its increase induces the well known exponential decrease of T'schroedLger' i n both 
the relativistic cases - for fixed U - the higher the value of L, the faster the oscillations of 
jrect bar ag a f unc ti on Q f £ Fig. s [5] and [6] illustrate this statement. 

To make our point we have now to take into account how a 'realistic' experiment would 
be carried out and to give the numerical values of the scales involved. 

Taking L ^> 1/m = Acompton = 0.024 A for electrons (and smaller for more massive 
bosons) is simply mandatory. In addition, in a 'real experiment', where almost monoener- 
getic particles - with average energy E and an energy uncertainty AE <C E - are gunned 
against the barrier, the in-state is obtained as some superposition of stationary states with 
E — \AE < E < E + \AE. Thus, even for reasonably monoenergetic electrons, a huge 
number of transmission resonances is involved in the diffusion process. As an example, for 
m = 0.51 MeV, U > 1.2 MeV, L > 10 A and E ~ m + \AE = 0.520 MeV (i.e. at the very 
beginning of the Klein zone with a relative monoenergeticity AE/E < 4%) no less than 
1450 spikes in the analogues of Fig.s[5]and[6]are involved. 
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Fi gure 5: ^ arr as a function of E in the Klein zone: U — 8m,, L — 4/m (up) and L — 6.5/?n 

(down) 




m Klein zone U-m 



1 


1 








m 


Klein zone 


U-m 



Figure 6: Xg 3Ctbarr as a function of E in the Klein zone: U = 8m, L = 4/m (up) and L = 6.5/m 
(down) 
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Such fast oscillations have to be taken into account [2]. This is done by replacing the 
sin 2 (. . . ) with | in (jSJ) and this, in turn, leads to an energy- averaged transmission coefficient 

T r ± ect bar = — L= (10) 

(the non-oscillating curve in Fig. [5] and El). As long as T ± is taken as an indicative 
prediction of the theory, the seeming paradox shows up in the following way: 



7=rect bar U~S>E I „ 

' S E 2 -m 2 2 f_ ( n ) 



E 2 -m 2 /2 m2 

indicating an 'unnaturally' high transmittivity of Dirac particles in presence of a high step. 
The 'unnatural' is referred to the comparison with the exponentially small Schrodinger 
prediction. Even without going to the first limit taken in ffTTj) . for example, for U = 2.5 
MeV and average kinetic energy E — m — 0.1 MeV, the 31.3% of the incident flux would be 
transmitted, the fraction going up to 64.7% for E — m = 0.5 MeV. 

The present paper addresses the question of whether Sauter's cure for the barrier is as 
good as it is for the step. In addition to giving the analytical result for 7^ auter barr ; our main 
result consists in showing that, for the Sauter barrier, in the part of the Klein zone where 
E — m ~ m, 

TpSauter barr ^ ^Sauter stcp^2 ^ e -2n m 2 £/U ^2) 

provided the slope £ is large enough so as to satisfy (JSJ). 

As a matter of fact, for devices involving either macro or mesoscopic dimensions, the 
result (TI21) is a doom: for practical purposes transmittivity is zero, in accordance with the 
naive expectation. Should one ever be able to set up a nanostructure with i as small as 
10 2_; " 3 Acompton and U not much higher than the threshold 2m, it is seen that the constraint 
dSJ) would be - even in this case - largely satisfied. Cases that could possibly be left uncovered 
by the present discussion are those of solid state physics where the 'effective mass' of the 
electron is smaller than the in- vacuum- value adopted here. In particular the case of graphene, 
where a vanishing effective mass is advocated by the authors of Ref. {f|, looks to date as 
the most promising ground where to observe the Klein paradox at work and, in the most 
optimistic case, to set up a "graphene device electronics" 

The paper is organised as follows. 



In Section [TT] the transfer matrix M step for a generic potential step is defined and the 
corresponding transmission coefficient discussed. In Section III the tranfer matrix and trans- 
mission coefficient for the corresponding space-inversion invariant barrier are expressed in 
terms of the matrix elements of M step . Section IV (relying on the analytical results extracted 
from Sauter's original paper and collected in the Appendix) illustrates the case of Sauter's 
trapezoidal barrier and justifies why the averanging procedure described for (fTDl) applies and 
how is it that (fT2l finally comes about. 



II. STEPS 



The present section contains material that is well known: as it mainly serves to introduce 
our notation, some of the statements to be found below will be made without proof. For 
the sake of conciseness we give a unified treatment of the Klein-Gordon equation in the 
first-order formalism and of the Dirac equation, in one space dimension. In the latter case 
we make use of two-component spinors since the spin - conserved by minimal coupling - is 
irrelevant. We will make use of the Pauli matrices: 

(o l\ (o -i 

\1 OJ \i 

as well as of the 2x2 identity matrix /. 
We are interested in the stationary states 

* E (x,t) = e- lEt ^ E (x), ^ E (x) = I rm I (13) 





of the Klein-Gordon one-dimensional wave equation 

tp' E = Xe 



(14) 



X ' E = -({E-U(x)) 2 -m 2 j M*) 
as well as of the Dirac equation. In the Pauli representation the one-dimensional Dirac 
Hamiltonian is 

= -1 a x 4- + 0-3 m + U(x) (15) 
ax 

i.e., spelling out the components, 

ip' E = i (E + m- U(x)) xe , 

(16) 

x' E = i(E -m-U (x)) ip E . 



U(x) = < 



(17) 



In both (Tl4|) and (|T6|) the energy E is understood to be in the Klein zone <^ and the step 
potential U(x) is given by 

x < 
W(a;) < a; < £ 

u x > e . 

The solutions of the above equations may all be expressed in the form 

au(p) e ipx + P u(—p) e~ ipx x<0 

E(x) < x < I 

ju(q) e tqx + 5u(—q) e~ tqx x > £ 

where S(rc), representing the solution in the intermediate region, evidently depends on the 
explicit form of U(x). The momenta p, q are given by (CQ), (EJ), and the two component wave 
functions are respectively 

/ 1 \ / 1 ^ 

, u D (±p) = 



= < 



(18) 



u K g(±p) 



\±<pj 
( i \ 



±p 
\E + ml 



(19) 



mkg(±?) 



\ ±lq J 



(20) 



( 1 

u D (±q) = 

M 

\E~U + mj 

Imposing the continuity of at the points x = 0, x = i entails two relations of linear 
dependence among the four constants a, (3, 7, 5. The transfer matrix M step , expressing 
such a dependence, is defined by 

7 ) = ( a ) = (" ") (" 

Use of the invariance of (}Tlj) and f|T6l) under charge conjugations 



(21) 



C: $(x) -> C**(x), C 



KG 



a 3 D 



(22) 



simplifies the form of M step down to 



M step 



a b 

b* a* 



(23) 



In addition, conservation of the currents (whose form is dictated by Noether's theorem), i.e. 
the Constance of 



with respect to x, implies 



where 



Jrp. 



|a| 2 -|/3p 



V± p 



KG 
D 

|7| 2 "|5| 2 ) 
+ = KG 



E-U + m 



E + m 

Finally, expressing 7, 5 in ff25l) as given by ff2Tl) . one obtains 



D . 



detMf p 



p -v ± 



(24) 



(25) 



(26) 



(27) 



In the Klein-Gordon case there is no doubt that the scattering state with the a source term 
to the left is identified by setting 5 = in (fTBj) . Indeed, according to (1251) . only an outgoing 
(from left to right) flux is left in the region x > i. The transmission coefficient is therefore, 
in general 

' " P 1 (28) 



rpStep _ 

KG — 



q 



v 



q \a\ 



Apparently Pauli was the first |l| to note the impact that P_ < has, in the Klein zone, on 
the identification of the scattering states. As the contributions to the flux of Dirac particles 
in the right end side of (125]) are interchanged with respect to the Klein- Gordon case, the 
transmission coefficient should be obtained by setting 7 = instead of S in (|18|) . This choice 
gives, for the generic form U{x) of the step, 



rpStep 



q E + m 



p X>_ 
q W 



p U — E — m 

to be contrasted with (I2"g|) . 

In the case of the rectangular step, explicit calculation yields 



(29) 



a 



rect step 
± 



2 V q 



ycect step 



5" 



(30) 



10 




Figure 7: j^ ectste P anc { j^tstep f or ^ _ x.2m as functions of U in the Klein zone 

whence 

± 4-X>± 

rprect step (? (31 1 

In the Dirac case it is evident that our identification of the scattering states is in disagree- 
ment, e.g., with the choice of Ref. [8(, where the interchange of the 7 and 5 terms is not 
effected. Our result is indeed obtained from theirs by effecting the substitution r — > — r (our 
p/q T>_ equals Bjorken-Drell's 1/r), which turns their transmission coefficient - negative in 
the Klein zone - into our positive expression. 

Considering now that, for fixed E, 



it follows that 



V ^ E — m 
lim -- £>_ = */— 32 

u^oo q V E + m v ' 



lim T* ect step = (33) 
u-,00 D E + p K J 



i.e. a finite value, to be contrasted with the Klein-Gordon case where 



entails 



lim V - = (34) 

u^oo q 



lim T^ tstcp = 0. (35) 

U— >oo 



The situation is summarized in Fig. [71 The comparison of the two cases is self evident. 
Whereas, in the first case, the barrier is transparent for essentially all the values of the 
potential, in the second one, apart from a region just above the threshold of the Klein 
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zone, the barrier becomes again impenetrable. This shows how wrong is the naive idea that 
only the oscillatory behaviour of the eigenfunction - due to the relativistic energy-momentum 
relation - does all the job in the Klein zone. Indeed the equation for the large component 
of the Dirac equation (obtained by taking the derivative of the first equation in (I16p and 
there using the second) is not equivalent to the Klein-Gordon equation (j!4p . The derivative, 
acting on U(x), yields extra terms e£ = —U'(x), thus supporting the expectation that a 
sharp edge is fundamental in making the difference between the previous results. This is 
indeed the content of Sauter's work [3j. 

At this stage we can no longer help interpreting the scattering states. Indeed, in the 
scattering state identified by choosing 7 = in (TT81) . the positive (i.e. toward the right) 
current in the region x > cannot come from electrons with negative kinetic energy, see 
([6]). It is rather due to positrons propagating backward in time, see viceversa (J7J). Hence - 
in the language of second quantization - a physical positron above the potential —U comes 
from the right, i.e the S term acts as a source of positrons available at x — +00. In this 
sense we are, therefore, in the presence of a e + -e~ annihilation process, and the reflection 
coefficient R = 1 — T necessarily is smaller than 1. Had we had alternatively chosen 5 = in 
(j!8p . we would have had, also from the right, an incoming current with an ensuing reflection 
coefficient bigger than 1. This would correspond to a flux of positrons travelling to the right, 
i.e. a process of pair creation. And this is the result much quoted as the 'Klein paradox', 
a statement that simply is not present in the original Klein's paper. Let us however again 
stress that even for R < 1 there is a paradox due to a finite fermion transmission through a 
high repulsive potential. Analogous considerations apply in the Klein-Gordon case with the 
previous provisions (Eq.s (|28|) and (!30|) and Fig. [7]) due to the reversed role of the 7 and S 
terms. 

III. BARRIERS 

The space-inversion invariant potential barrier obtained by using twice the step (ITTj) is 

W(x) = U(x + £ + L) +U(-(x + £ + L)) (36) 
and the corresponding stationary states have the form 
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au{p)e ipx + f3u{-p)e~ ipx x < -{£ + L) 

S_(z) -(£ + L) < x < -L 

-iu(q)e iqx + 5u(-q)e- iqx -L < x < L (37) 

E + (x) L < x < £ + L 

a'u(p)e ipx + p'u(-p)e~ ipx x>£ + L 

(the coefficients a, /3, 7, 5 appearing in the latter equation should not be confused with those 
appearing in f|T8l) to which they are however related through the translation x — > x + i + L) . 
Again, the continuity relations at x = — (I + L), —L and x = L, £ + L respectively entail 



7 



M 



left 



a 



P. 



a' 
3' 



M 



right 



7 



(38) 



whence the transfer matrix relative to the barrier 



a \ , 1 at 

= M bar 



.P. 



P. 



A B 

B* A* 



a 



is obtained as 



M bar = M right M left . 



(39) 



(40) 



The connection of M right and M left with M step is established by taking translations and 
reflection into account. The effect of the translation to the left x — > x + £ + L is expressed 
with the aid of the matrix 

( e ika q 
e~ lka 

One gets 



(41) 



M 



left 



T(q,£ + L) M step T(-p,£ + L) 



(42) 



whereas M right is obtained either by reflecting x — > —x in M left , or by first reflecting M step 
and then translating to the right x — > x — (£ + L): 



M right = ^ ( M Mtyl ai = 

T{p, -(£ + L)) (7 1 (M ste P)-V 1 T( - q, -{£ + L)) 



(43) 
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Exploiting (|40|) - (H3|) one finally obtains for the coefficients appearing in (1391) 

a 2 e 2i(q-p)(e+L) _ b *2 e -2i(q+p)(e+L) 

\a\ 2 ~\b\ 



u c - - ' — u e ^ - - 

= r^-n^ ( 44 ) 



B = («) 

In the case of the barrier we do not have the ambiguity - connected with the identification 
of the scattering states - we have discussed for the step in ( 1281) and ( 1291) . The transmission 
coefficient of the barrier is 1 / 1 ^4 1 2 both in the Klein-Gordon and the Dirac case. When it is 
expressed in terms of the two matrix elements entering the transfer matrix of the step 

a=\a\e^, b=\b\e ir >, (46) 



takes the form 



Tbar = 4HW " (47) 

1 + (| |2-| 6 |2)2 81112 ^ (t + L)+£ + v) 



This can be further simplified on account of (l2Tl) - (l29l) : 



-ibar 



ystep 

1 + 4^ F ^sin 2 (2 q (£ + L)+^ + r ] ) 



(48) 



The above formula is consistent with the results cited in the introduction, i.e. when (|28|) - (|3ip 
- appropriate for the rectangular barriers - are used, ([8]), (Q are reobtained. 

The dependence of T bar on the energy is through the phases £, rj of the transfer matrix 
elements f )46|) . as well as through T step itself and q. The dependence on the barrier width 
2L is, instead, completely spelled out in the argument of the sin function. This means that 
for given energy, no matter how opaque the step U(x) may be, the width parameter L can 
be always adjusted in such a way that the corresponding even barrier is totally transparent 
(this is analogue to the situation of dielectric stratified media in optics: see e.g. js], where 
however only the case of sharp edge I = 0, i.e. £ = rj = 0, is considered). Going back to 
particles, as discussed in the introduction, there may be the difficulty that, for given average 
energy E and energy uncertainty AE of the incident particles, any 'realistic' value of L 
may be such that too many maxima of T bar , due to the oscillations of sin 2 (. . . ), come into 
the range (E — ^AE,E + ~AE). The averaging procedure sin 2 (. . .) — > | in fl4"8"j) is then 
necessary and f|T2|) follows immediately when T step 1. 
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Figure 8: xg autbar for U = 3 m, L = 0, £ = 3 A C om P ton as function of E in the Klein zone 
IV. SAUTER TRAPEZOIDAL BARRIER FOR THE DIRAC CASE 

The analytic expression for ySauter bar can ^ e (j e( j uce( j f r0 m the long formulae we have 
extracted from Sauter's paper and reported in the Appendix. The graphs below are more 
eloquent than the expression. In assigning numerical values to the parameters, we have in 
mind the case of electrons in vacuum. Therefore the unit of energy is m = .51 MeV and the 
unit of lenght is Acompton = -024 A. 

We repeat here the warning already made in the Introduction: there may be cases where 
the units change by orders of magnitude and the following discussion may not apply. 

In consequence of the scales we have chosen, the values we give to U are not much higher 
than the threshold, typically U = 3 m. As for i and L, certainly they both have to take a 
value of the order of at least 10 2 Acompton- It is however instructive to vary one parameter 
at a time. 

Fig. [S] plots J , p autbar against the energy E in the Klein zone for U = 3m, £ = 3Ac mpton 
and L = 0: it is a triangular barrier that shows the effect of I alone. It indicates that a 
flux of monoenergetic electrons (AE <ti width of the peak) is entirely transmitted. This 
may come a little bit as a surprise because, referring to Fig. [2] and in a pictorial space-time 
description, the incoming electron encounters: 

(i) a free propagation region (from x = — oo to x — 0); 

(ii) a region from to x\ that is still classically accessible, but where the propagation 
is non-free: the two solutions behave as (think e.g. of the semiclassical approximation) 
exponentials of opposite imaginary arguments; 



15 



(iii) the classically inaccessible region (xi,^) with non-propagating solutions: the two 
solutions of the Dirac equation behave as exponentials of opposite real arguments (a region 
that we will improperly call the 'attenuation' region); 

(iv) eventually the classically inaccessible region from x 2 to I where propagation, albeit 
non-free, still occurrs; 

then, from t to +00, the same steps just described, but in reverse order. 

If one takes the Schrodinger case as a guidance, one might expect attenuation of the wave 
function both in the left and in the right attenuation region. This happens for almost all the 
values of E in the Klein- zone, with the exception of the peak displayed. Indeed, for average 
energy E at the center of the peak, the single monocromatic components, that make up the 
incident packet, undergo two transmissions at the edges and possible multiple reflections in 
between, but they succeed in keeping somehow memory of the relative phases and the packet 
is (almost entirely) reconstructed beyond the barrier. For the values we have indicated the 
'device' could serve as a monocromator for incident wave packets with energy uncertainty 
greater than the width of the displayed peak. From a practical point of view this is illusory, 
due to the exceedingly small value of I = 3Acom P ton = 0.072 A. 

Fig. [9] shows the impact of increasing L from to lOOAcompton = 2.4 Ain the preceding 
case. A region of strict free propagation is interposed, the peak of Fig. [8] is replicated several 
times (the argument of the sin 2 (. . . ) in (T48]) oscillates much faster) and the peaks displayed 
all reach the value 1. However we have chosen to reduce the scale of the ordinates to make 
the energy averaged y^ autbar visible. The oscillations of the non-averaged T are so fast (and 
still the value given to L may be largely considered too small) that only T is susceptible of 
comparison with the data of whatever 'Gedanken' experiment. Nonetheless a value of T of 
the order of some part per thousand could still be considered an interesting signal. 

Once the role of L is understood, in Fig. [10] we go back to the case L = but raise 
the value of £ from 3 Acompton to 100 Acompton- The main difference is in the scale of the 
ordinates. Again the peaks all reach the value 1, but the scale of T has dropped to io~ 90 ^ 91 . 
As for the part of the Klein zone that also satisfies E — m ~ m, this is in perfect agreement 
with (jT2j)). 

The analytical result, displayed in the figure, shows that the estimate extends to the 
whole Klein zone. 

The number of peaks in the Klein zone increases with increasing £, but the width of each 
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Figure 9: Jg autbar for U = 3 m, L = 100 A C om P ton, I = 3Ac m P ton as function of E in the Klein 
zone 




m=l U-m=2 



Figure 10: T^ autbar for U = 3 m, L = 0, i = lOOAcompton as function of E in the Klein zone 

of them shrinks much faster: this is why the value of T has dropped of so many orders of 
magnitude. 

Turning L on will not appreciably change the scale of the ordinates, but the energy 
averaging becomes unavoidable. This washes out any phase information among the mo- 
noenergetic components the incident packet is made of. As a consequence the packet, that 
has undergone an attenuation in the attenuation region at the first edge, may only undergo 
a second attenuation at the second edge. The result is that the transmission coefficient (|48|) 
is shattered down to the asymptotic value given in ffl2|) . 

Finally the trapezoidal form of the barrier, we have chosen in order to be able to exhibit 
analytical results, should not be crucial to the above conclusion. The same qualitative 
result should obtain with all the barriers where the dimension I of the regions that support 
a relevant electric field (eS > 2m 2 ) is large enough so as to fulfill <^j. Ref. [nj], where again 
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Sauter provides analytical results for the potential 



as well as the relevant asymptotics, substantiates the above expectation. 

The foregoing ar gum ent should also work for the 'derealization' of particles in supercrit- 



ical potentials [111 ] . [12] . This will be possibly considered elsewhere. 
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Appendix: LINEARLY RAISING POTENTIAL AND SAUTER SOLUTION FOR 
THE DIRAC CASE 

Referring to the Dirac equation with the step potential of Fig. [2] 

U{x) = Uj, 0<x<£ (A.l) 

the solution in the intermediate region < x < £ to be fed in 0181) can be worked 
out from Sauter's article [3(. His representation for the Dirac Hamiltonian is different from 
Pauli's f|T5l) : he uses 

^Sauter = { + N _ a ^ m = y #Pauli y\ 

dx 

with (up to an irrelevant phase factor) 

V = i(J-i>i + <7 2 -<7 3 )) • (A.3) 



He then introduces the argument 



&) = \T(*-7t*) ( A - 4 ) 



u 



and the degenerate hypergeometric function 



Ha ^ z) = l + -- + J (f -1- + ... (A.5) 
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in terms of which the two independent eigenfunctions belonging to energy E 

X; • | | X] • CE = • j ' } (A.(») 

k 9e{x) 

in the intermediate region < x < I are given by 

,,(,)=<<«■/» •gl2L£, (A.7) 

9E (x) = -m Jle'C* 4> (1 ^ + 1, \, -i f ) . (A.8) 
Thanks to (1A.3I) . the to be fed in ffl8l) is the linear combination 

, Mx)+i^(x)\ +i/ / 9%(x)+if* E (x) 

The transfer matrix M D aut step is obtained by imposing the continuity of ([18]) at x = and 
x = i and then eliminating /i and v. Its matrix elements (123]) turn out to be 



p-iqt 

^Saut step 



Aq (E + m) 

iql 

Aq (E + m) 



(p+Q+R - P-Q-R* + % {P+QS + P-Q+S*)) (A.10) 



6 Saut step = (p_Q + R - P + Q_R* + i (P_Q_S + P+Q+S*)) (A.ll) 

Aa ( E + m ) V / 
where 

R = f E (0)fm-9EiO)9E^) (A.12) 

5 = / e (0)^)-5 E (0)/eW (A.13) 

P± = E + m±p (A.14) 

Q± = E-U + m±q. (A.15) 

Substitution of the above formulae in (|4TJl) - (j4"81) provides the analytic expression for T^ aut stop 
by which the graphs of Fig.s IHlfTOl above are obtained. Finally Sauter himself works out the 
asymptotic form (TJJ that T^ aut step has when both E — m ~ m and are satisfied: 

T Saut step = q(U -E + q) e ^ mH/U , x 0{U/m 2^ 

having included a prefactor omitted in ([T2]) . 
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